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It is shown that the DAMA data indicate two dark matter components, one that circulates around 
the galactic center (GC) and another that is emitted from the GC. From the location of the maximum 
yearly variation, one can compute the ratio of the two components. 
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I. INTRODUCTION 

The yearly variation observed in the DAMA 
experiment [H,Q extends for 12 years and the location 
of the maximum has been established with 8.9a. Since 
the time at the maximum does not coincide with that 
when the Earth orbit approaches the GC or is tangential 
to a circular orbit around the GC, it is natural to assume 
that dark matter has two components, one that circu- 
lates around the GC and another that is emitted from 
the GC. This article computes the ratio of the two cm- 
ponents, and explores its significance. A possibility of the 
Migdal effect being a source that differentiates the results 
of the DAMA data and the CDMS data is discussed. 



II. TWO DARK MATTER PARTICLE (DMP) 
COMPONENTS 



If the GC is on the 18h of the RA coordinate, the 
Earth is moving towards it on March 22nd (81 days from 
January 1st) and moves tangentially to a circular orbit 
around the GC on June 22nd (173 days) and December 
22. The real coordinates of the GC areQ 

RA = 17 /l 45 m 37M991 = 266°.40499625 (1) 
DEC = -28°56'10".221 = -28°.93617242 (J200©,) 

then the time when the Earth is moving towards the GC 
is March 18.35 or 

<oi = 77.35 in days (3) 

and the time when the Earth is closest to the GC is 

f,Q2 = 168.66 in days. (4) 

Since neither of these dates coincides with the observed 
time of maximum yearly variation [l|, 



to = 145 ±7, (8.9o-), 



(5) 



one has to assume that there are two dark matter com- 
ponents. Let us denote the yearly intensity variation of 
the component coming from the GC by A\ and that of 



the component along a circular orbit by A 2 . The total 
DMP intensity variation that yields A is given by 

Ao,o%(u(t — 1 )) = Ai cos(w(i — ioi)) + A 2 cos(uj(t — to2)), 

(6) 

which gives 



Acos(u;to) — Ai cos(ujt 01 ) + A 2 cos(ujt 02 ) (7) 



and 



where 



Asin(wto) = A\ sin(o;ioi) + A 2 sin(o;io2), 



2tt 



{in days ) 



(8) 



(9) 



365.2422 

From Eq. (0 and Eq. ((8]), one gets 

A 2 = A\+Al + 2A 1 A 2 cos(u{t m -t 01 )) (10) 

(11) 



A( + Ai, 



since 



^(^02 — *oi) = 2> 



and 



Ai sin(^t i) + A 



- 2 sin (u)t 02 ) 

Ai cos(wioi) + A 2 cos(tut Q2 ) ' 

Solving Eq. (Tl3|) for A 2 /Ai, one obtains 

A 2 _ sin (a; fa - t m )) 
Ai sm(u(t Q2 - t )) ' 



(12) 



(13) 



(14) 



2.35 ± 1.04. 



(15) 



From Eq. (|3][5j> , one gets 

M 

Ax 

The DMP component that circulates on a circular orbit 
around the GC is moving with the solar system and it 
creates a yearly variation due to a variation in the radial 
DMP distribution, as is described in the next section. 
This gives the A 2 component. In March, when the Earth 
is approaching the GC, it encounters DMP emitted from 
black holes at the GC with higher relative velocity than 
in September, when the earth is moving away from the 
GC, with smaller relative velocity. This gives the A\ 
component. The details of the both components will be 
described in the following sections. 
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III. DMP CIRCULATING THE GC 

The simplest assumption that gives the flat velocity 
distribution observed in most spiral galaxies is a spheri- 
cally symmetric distribution of DMP. If one assumes that 
DMP are orbiting in all possible directions, no yearly 
variation is produced, since a summer-winter difference 
is cancelled by two DMP orbitting in the opposite direc- 
tions. A yearly variation is produced by a radial distri- 
bution of DMP. If the DMP distribution is a decreasing 
function of distance from the GC, then the Earth en- 
counters more DMP flux in the summer due to the fact 
that Earth is closer to the GC in the summer. This is 
consistent with the DAMA data in the previous section. 

If the orbital motion of DMP is constrained so that 
the angular momentum component relative to that of 
the solar system or the visible spiral galactic plane, then 
the component of the orbital motion in the galactic plane 
contributes to the yearly varition, since the component 
vertical to the galactic plane cancels between orbits with 
opposite vertical components. Assuming that the average 
orbit is circular, its component in the galactic plane has 
a velocity less than that of the solar system. Then the 
maximum of yearly variation occurs in the winter, since 
the orbital motion in the vertical plane is opposite to that 
of the solar system relative to the latter, in contradiction 
to the computation in the previous section. In order to 
get a solution with a maximum in the summer, one has to 
have a contribution from a decreasing radial distribution 
compensating the effect of orbital motion. 

In both cases, the effect of the radial distribution has 
to dominate the yearly variation. 



IV. EMISSION OF DMP FROM GC BLACK 
HOLES 

The existence of DMP in orbital motion around the 
GC is quite natural. It is also natural to assume that 
its distribution is a decreasing function of distance from 
the GC at the location of the solar system, in order to 
get a flat velocity distribution around the GC. However, 
the existence of DMP emitted from the GC requires an 
explanation. That is related to a recent discovery by the 
Pierre Auger Project concerning the correlation between 
high energy cosmic rays and the location of AGN (Ac- 
tive Galactic Nuclei) [J] and a model of the author which 
predicted the Auger data since 1985. 

The Pierre Auger Project data suggest that high en- 
ergy cosmic rays are emitted from AGN, massive black 
holes. In a series of articles (5I] - fl3j ■ the author has pre- 
sented a model for the emission of high energy particles 
from AGN. The following is a summary of the model. 

1) Quantum effects on gravity yield repulsive forces at 
short distances 0|, 0. 

2) The collapse of black holes results in explosive 
bounce back motion with the emission of high energy 
particles. 



3) Consideration of the Penrose diagram eliminates the 
horizon problem for black holes [8j. Black holes are not 
black anymore. 

4) The knee energy for high energy cosmic rays can 
be understood as a split between a radiation-dominated 
region and a matter-dominated region, not unlike that in 
the expansion of the universe. (See page 10 of the lecture 
notes ll-QO 

5) Neutrinos and gamma rays as well as cosmic rays 
should have the same spectral index for each AGN. They 
should show a knee energy phenomenon, a break in the 
energy spectral index, similar to that for the cosmic ray 
energy spectrum. 

6) The recent announcement by Hawking rescinding an 
earlier claim about the information paradox[l4j is consis- 
tent with this model. 

Further discussion of the knee energy in the model pre- 
dicts the existence of a new mass scale in the knee energy 
range, in order to have the knee energy phenomenon in 
cosmic rays [15j . The following are additional features of 
the model. 

7) If the proposed new particle with mass in the knee 
energy range (0.1 PeV~2 PeV) is stable and weakly in- 
teracting with ordinary particles, then it becomes a can- 
didate for a DMP. It does not necessarily have to be a 
supersymmetric particle. That is an open question. How- 
ever, if it is supersymmetric, then it is easy to make a 
model for a weakly interacting DMP[16]. The only re- 
quirement is that such particles must be present in AGN 
or black holes so that the the knee energy is observed 
when cosmic rays are emitted from AGN. A suggested 
name for the particle is sion (xion) (l5j . using the Chi- 
nese/Japanese word for knee, si (xi). 

8) If the particle is weakly interacting, then it does 
not obey the GZK cutoff [13], since its interaction with 
photons in cosmic backgroud radiation is weak, as was 
pointed out earlier. This is a possibe resolution of the 
GZK puzzle[I3, 03, & 

In summary, this model predicted the Pierre Auger 
Project data. Moreover it suggests the existence of a 
new particle in the PeV mass range, in order to explain 
the knee energy phenomenon in the cosmic ray spectrum. 
The author has suggested that this new particle is a can- 
didate for DMPglj. 

V. AGN AND GC EMITS DMP 

We assume that the incident particles above the 
GZK cutoff observed by the Akeno-AGASA detector are 
weakly interacting particles at the PeV mass scale, which 
are required to exist in order to explain the cosmic ray 
knee energy in the model. In a model where the accel- 
eration takes place by gravity such as that proposed by 
the author, there is no difficulty in accelerating a weakly 
interacting and neutral DMP by gravity. 

Then AGASA and DAMA data present a consistent 
picture for DMP emission from black holes. The only dif- 
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ference is the AGASA data captures the high energy end 
of the DMP spectrum, while the DAMA data captures 
the lower end: In the latter, the cross section increases 
linearly with energy and the energy spectrum of the par- 
ticle distribution decreases with energy, E~ 2 5 ~ E~ 3 , so 
that the lower end of the spectrum contributes to the 
DAMA data. 

The discussion in the previous section and this sec- 
tion provides a basis for the DMP component emitted 
from the GC that is implied by the DAMA data. Due 
to the nature of weak interactions, where cross sections 
are an increasing function of energy, March has a higher 
probability of encountering with the dark matter than 
September. 

VI. THE MIGDAL EFFECT AND THE DAMA 
DATA 

An important question is why other experiments such 
as CDMS have not observed a yearly variation of DMP. 



It was suggested that enhancement by the Migdal effect 
might be a key factor f-^, @. In ref. [12], the Migdal 
effect is computed for a Nal detector and it is shown that 
there is a large ionization effect for the inner electron. 
The Migdal effect is known to give a smaller effect for 
targets with large atomic number Z, since the Coulomb 
binding is proportional to (Ze) 2 and therefore the Migdal 
effect is proportional [23| to 1/Z 2 . The Migdal effect on 
the Ge target of CDMS is 1/ 10 of that of Na target of 
DAMA. That might explain the different results for the 
two detectors. 
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